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Abstract. For the spatially homogeneous Boltzmann equation with cutoff hard 
potentials it is shown that solutions remain bounded from above, uniformly in 
time, by a Maxwellian distribution, provided the initial data have a Maxwellian 
upper bound. The main technique is based on a comparison principle that uses 
a certain dissipative property of the linear Boltzmann equation. Implications 
of the technique to propagation of upper Maxwellian bounds in the spatially- 
inhomogeneous case are discussed. 

1. Introduction and main result 

The nonlinear Boltzmann equation is a classical model for a gas at low or moderate 
densities. The gas in a spatial domain Q C R d , d > 2, is modeled by the mass density 
function f(x,v,t), (x,v) G f2 x R , where v is the velocity variable, and t G R is 
time. The equation for / reads 

(1) (dt + v-V x )f = Q(f), 

where Q{f) is a quadratic integral operator, expressing the change of / due to instan- 
taneous binary collisions of particles. The precise form of Q(f) will be introduced 
below, cf. also [10,34]. 

Although some of our results deal with more general situations, we will be mostly 
concerned with a special class of solutions that are independent of the spatial variable 
(spatially homogeneous solutions). In this case / = f(v,t) and one can study the 
initial-value problem 

(2) d t f = Q(f), f\ t=0 = fo, 

where < /o G L 1 (JBL d ). The spatially homogeneous theory is very well developed 
although not complete. In the present paper we shall solve one of the most noticeable 
open problems remaining in the field, by establishing the following result. 

Theorem 1. Assume that < fo(v) < M (v), for a. a. v G M. d , where M (v) = 
e -aoM +c £ S ^ e d ens ity j a Maxwellian distribution, oq > 0, cq G R. Let f(v,t), 
v G M. d , t > be the unique solution of equation (J2]) for hard potentials with the 
angular cutoff assumptions ([5]), ([7]) ; that preserves the initial mass and energy ( 1121) . 
Then there are constants a > and c G R such that f(v, t) < M(v), for a. a. v G R d 
and for all t > 0, where M(v) = e - a \ v \ 2+c . 

Before going on, let us make a few comments about the interest of these bounds. 
Maxwellian functions 

M(v) = e - a H 2 +^+ c ; W ith a > 0, c G R, b G R d constants, 
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are unique, within integrable functions, equilibrium solutions of (j2J), and they pro- 
vide global attractors for the time-evolution described by (J2J) (or (jTJ), with appro- 
priate boundary conditions). Classes of functions bounded above by Maxwellians 
provide a convenient analytical framework for the local existence theory of strong 
solutions for ([!]), see Grad [22] and Kaniel-Shinbrot [25]. Such bounds also play an 
important role in the proof of validation of the Boltzmann equation by Lanford [27], 
see also [10]. However, establishing the propagation of uniform bounds is generally 
a difficult problem, solved only in the context of small solutions in an unbounded 
space, see Illner-Shinbrot [24] and subsequent works [4,21,23,29]. The above results 
rely in a crucial way on the decay of solutions for large \x\ and on the dispersive effect 
of the transport term, in order to control the nonlinearity. These effects may not 
be significant in other physical situations, and the spatially homogeneous problem 
presents a simplest example of such regime. 

In the spatially homogeneous case many additional properties of solutions can 
be established. Upper bounds with polynomial decay for \v\ large hold uniformly 
in time, see Carleman [8,9] and Arkeryd [2]. Solutions are also known to have a 
lower Maxwellian bound for all positive times, even for compactly supported initial 
data [32]. Many results have been established that concern the behavior of the 
moments with respect to the velocity variable, following the work by Povzner [31], 
see in particular [1,6,12,15,30]. The Carleman-type estimates [2,8,9] were crucial 
in the treatment of the weakly inhomogeneous problem given in [3]. However, as 
also pointed out in ref. [3], Maxwellian bounds of the local existence theory [22,25] 
are not known to hold on longer time-intervals, and it remains an open problem 
to characterize the approach to the Maxwellian equilibrium in classes of functions 
with exponential decay. The present work aims to at least partially remedy this 
situation, and to develop a technique that could be used to obtain further results in 
this direction. 

We will next introduce the notation and the necessary concepts to make the 
statement of Theorem [1] more precise. We set in (j2j) 



(3) Q(f)(v,t)= / (flf-f*f)B(v-v*,a)dadv*, 

JR d J 5 d -! 

where, adopting common shorthand notations, / = f(v,t), f = f(v',t), /* = 
f(v*,t), fl = f(v'„t). Here v, v* denote the velocities of two particles either before 
or after a collision, 

,. s , v + v* \v — vJ . v + v* \v — vJ 

4 v — 1 a, t> = a, 

\ ) 2 2*2 2 

are the transformed velocities, and o G S d ~ x is a parameter determining the direction 
of the relative velocity v' — v[. In the more general case of ([T]), the space variable x 
appears (similarly to t above) in each occurrence of /, /#, /', /£; we shall often omit 
the t and x variables from the notation for brevity. 

Many properties of the solutions of the Boltzmann equation depend crucially 
on certain features of the kernel B in (J3j). Its physical meaning is the product 
of the magnitude of the relative velocity by the effective scattering cross-section 
(see [26, §18] for terminology and explicit examples); this quantity characterizes 
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the relative frequency of collisions between particles. Our assumptions on B fall in 
the category of "hard potentials with angular cutoff", cf. [34]. More precisely, we 
assume that 

(5) B(v — v m , a) = \v — v*\P h(cos$), cos^ = ^~"*' > | Cr , 

where < (3 < 1 is a constant and h is a nonnegative function on (—1,1) such that 

(6) h(z) + h(—z) is nondecreasing on (0, 1) 
and 

(7) < /i(costf)sin a tf < C, ■& E (0,tt), 

where a < d — 1 and C is a constant. Assumption (J7J) implies in particular that the 
integral j sd -i ftr(cosi?) da is finite; for convenience we normalize it by setting 

(8) / h(cos$)da = uu d „ 2 [ h{z) (1 - z 2 )^ dz = 1, 
Js*- 1 J-i 

where uJd-2 is the measure of the (d — 2)-dimensional sphere. The classical hard- 
sphere model in M. d , satisfies (jSJ) with (3 = 1, (jSJ) and (j7j) with a — d — 3. 

Notice that we can write Q(/) = Q + (f) — Q~(f), where Q + (f) is the "gain" term, 
and Q~(f) is the "loss" term, 

Q + (f)= [ I f'f:B(v-v*,a)dadv*, Q~ (f) = (/ * \vf) f, 

JR d JS 11 - 1 

and * denotes the convolution in v. Because of the symmetry o \— > — a in the integral 
defining Q + {f) we can restrict the a- integration above to the half-sphere {cos?9 > 0} 
if we simultaneously replace B{v — u*, a) by 

S(u - v*, a) := (B(v - u*, a) + S(u - u*, -cr)) l{ cos tf>o}- 

It will be convenient to introduce the following (nonsymmetric) bilinear forms of the 
collision terms, 

(9) Q + (f,g)= [ [ f:g'B(v-v*,a)dadv^ Q~ tf , g) = tf * \vf) g , 

for which obviously Q ± (f) = <5 ± (/, /)• 

We say that a nonnegative function / £ C([0, oo); L 1 (R d )), such that (1 + |f | 2 )/ G 
L°°((0, oo); L 1 (R d )), is a (mild) solution of (0) if for almost all v eR d 

(10) /(«,0) = / (v); f(v,t)-f(v,s)= I Q(f)(v,r)dr, 

J s 

for all < s < t. Notice that the conditions on / imply (in the spatially-homogeneous 
case!) that 

(11) g + (/),g-(/)GL oo ((0,oo);L 1 (R d )), 

so the integral form in ( TlDl is well-defined. This also implies that / is weakly 
differentiable with respect to t and that the differential equation (T5]) holds in the 
sense of distributions on R d x (0, oo). 
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The existence of a unique solution satisfying the conservations of mass and energy, 
(12) / f(v,t)dv= [ fo(v)dv, [ f(v,t)\v\ 2 dv= [ f (v)\v\ 2 dv 



follows from a theorem by Mischler and Wennberg [30], for all fo > for which 
the above integrals are finite. The second condition in ( fl2l) is also necessary for the 
uniqueness [35] . For the initial data with strong decay (as in Theorem [1]) one could 
also refer to the well-known results by Carleman, Arkeryd and DiBlasio [1,2,13]. 

The following theorem summarizes the main results about qualitative properties 
of solutions in the case of "hard potentials with cutoff" known before this work. 

Theorem 2. Let f(v,t), v £ R d , t > 0, (n > 2) be a solution of (|2J) that satis- 
fies (I12p . and let the kernel B in the Boltzmann operator ([3]) satisfy (151) . (I7j). Then 

(i) if f £ L°°(R d ) then /(t, •) £ L°°(R d ), t > 0. Moreover, if (1 + |u|)'/o £ 
L°°(R d ) for some s > s , then (1 + \v\) s f(v, t) £ L°°(M. d ), t > 0. Here sq is a 
constant dependent on the dimension d. 

(ii) if the integral of f is nonzero, then for every to > there is a Maxwellian 
M(v) = Ke-^ 2 , K > 0, k > such that 

f(v, t) > M(v), t > t , for a.a.ve R d . 

(iii) for all t > and for all k > 1, the quantity rrik{t) = j Rd f(v,t) \v\ 2k dv 
is bounded uniformly for t > to; moreover, this bound is uniform in t > if 
m fc (0) < +oo. 



(iv) Inthe case d = 3 and B(v — v*,cr) = c\v — v*\ (hard spheres) or B{y — v*, a) = 
h( h £ L 1 (— 1, 1) (pseudo-Maxwell particles) if fo satisfies 

A e LV) 

for some Maxwellian Mq(v) = e~ a °^ 2 , ao > 0, then there exists constants a > 0, 
C such that 

f(v,t) 



M(v) 



dv < C, 



where M(v) = e a \ v \ 2 . 



Part (i) of this theorem is due to Carleman [9] in the case of the hard spheres; 
the general case was studied by Arkeryd in [2]. Part (ii) is due to A. Pulvirenti and 
Wennberg [32]. Part (iii) is due to Desvillettes [12] under the additional assumption 
that a moment rrik (t) of order k > 1 is finite initially; this assumption was removed 
by Mischler and Wennberg [30]. Earlier result by Arkeryd [1] and Elmroth [15] state 
that all moments remain bounded uniformly in time, once they are finite initially. 
Finally, part (iv) is due to Bobylev [6]; we will give an extension of this result to 
the class of Boltzmann kernels satisfying (151)— (17)) in Section [2j 

Our main contribution in the present work is to show that the estimates for the 
spatially homogeneous Boltzmann equation (precisely, parts (i) and (iv) of Theo- 
rem [2], together with the conservation of mass) imply Theorem [TJ Since we do not 
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use other properties of the spatially-homogeneous problem we can state our result 
in a more general, spatially inhomogeneous setting. 

We consider solutions of flTJ with the spatial domain Q = T d (rf-dimensional torus, 
or the unit hypercube with periodic boundary conditions), on an arbitrary finite time 
interval [0, T\. Spatially homogeneous solutions are then a special subclass charac- 
terized by the constant dependence on the x variable. To simplify the presentation, 
let us assume sufficient regularity (smoothness) of the solutions f(x,v,t) with re- 
spect to the x and t variables; this is not a restriction in the setting of Theorem 1, 
and the requirements of smoothness will be relaxed significantly later on to include 
a sufficiently wide class of weak solutions of the spatially inhomogeneous problem. 

Theorem 3. LetT > and let f £ C([0, T]; L l {T d x R d )), f > 0, be a (sufficiently 
regular) solution of the Boltzmann equation (CQ), with the initial condition 

f{x,v,0) = f (x,v) < M (v), for a. a. (x,v) eT d x R d , 

where M (v) = e~ a °^ +c ° , a > 0, cq E IBL Assume that the solution f(x,v,t) 
satisfies the estimates 



(13) / f(x,v,t)dv > po, (x, t) eT x [0, T], 
and 

(14) sup \\f(x,v,t)\\ L ~ < Co, sup / f^A dv<C Xi 

(x,t)eT d x[0,T] (x,t)£T d x[0,T] JR d M l{ v ) 

where M\(y) = e ~ ai l 1 'l 2+Cl and < a x < a 0; c 1; p , C , C\ are constants. Then for 
any < a < a\, for any t G [0, T] 

f{x,v,t) < M(v), for a. a. (x,v) e T d x R d , 

where M{v) = e - a \ v \ 2+c j and the constant c depends on a, a , c , a\, c\, p , C and 
C\ only. 

Remark. The regularity assumptions in Theorem [3] are not particularly restrictive. 
The precise conditions in the spatially inhomogeneous case are that / is a mild 
(renormalized) solution of (pQ) that is dissipative in the sense of P.-L. Lions (see 
Definition [TOl in Section [3]). A sufficient condition that is naturally satisfied in the 
spatially-homogeneous case is that (fTTI) holds in addition to ([TO]) . 



The plan of the paper is as follows. In Section [2] we extend property (iv) from The- 
orem [2] to the class of Boltzmann kernels satisfying (jSD^dZ])- This involves a rather 
technical analysis, specific to the spatially-homogeneous problem, which is based 
mainly on a development of the ideas from [5-7] . The result of Section [2j however, 
illustrates an important point that the type of behavior described by Theorem [1] is 
not a particular feature of the hard-sphere model, but rather a generic phenomenon 
that holds for a wide class of collision kernels of "hard" type. The key step occurs in 
Section [3j there we introduce the technique based on a comparison principle which 
plays a crucial role in the derivation of pointwise estimates. In Section H] we prove a 
weighted bound for the collision term, based on the Carleman representation of the 
gain operator, which is used in the comparison argument. Finally, some classical 
results used throughout the text are recalled in three Appendices. 
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Convention: Throughout the text, the function sign z is defined as 1 for z > 0, — 1 
for z < and an arbitrary fixed value in [—1, 1] for z — 0. 

2. Weighted L 1 estimates of solutions 

The aim of this section is the following result, originally due to Bobylev in the 
case of the "hard spheres" and Maxwell molecules [5,6]. 

Theorem 4. Let f(v,t), v G Mr, t > (n > 2) be a solution of the spatially 
homogeneous Boltzmann equation (j2J) with the collision kernel B satisfying ([5])-(j7]) 
and with the initial datum fo>0 such that 

(15) A G L\R d ) 
K ' M 

for a certain Maxwellian Mq(v) = e~ a °^ 2 , where ao is a positive constant. Then 
there exist constants C , a > 0, such that 

(16) / ^4dv<C, t>0, 



M(v) 

where M(v) = e~ a ^ 2 . 

Our approach to the problem is based on the analysis of the sequence of moments, 
(17) m k {t)= [ f(v,t)\v\ 2k dv, A; = 0,1..., 



and particularly, of the growth of mj.(t) as k — > oo. The relation between the 
moments (IT7|) and the weighted averages ffl6|) is given by the formal expansion 



(18) / ^4 dv = V a* 

In view of (fl8l) , to prove Theorem H] it suffices to show 

(19) lim sup fc - r = 0, for some A large enough. 

Our proof of (fl9l) is to a large extent a refinement of the original approach in [6] . 
One particular technical aspect which allows us to simplify some of the arguments 
is the systematic use of the interpolation inequalities 



(20) l^l V 1 < p^i y < (M kl <k<k 2 



m kl (t) \ H < f m k (t) \ k ^ f m k2 (t) \ k 2 
m o J ~ V m o 

which follow directly from (TTTf) by application of either Holder or Jensen's inequali- 
ties. 

It is well-known that if the kernel B(\v — t>*|, cosi?) in is constant in the first 
argument (the case of the Maxwell, or pseudo-Maxwell, particles) then the equations 
for the moments nikit) with integer k form a closed infinite system of ODE. This 
property no longer holds if the kernel B depends on \v — vJ, and one has to work 
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with inequalities instead of equations. If the kernel B has the homogeneity \v — vJ 13 , 
one also generally has to consider the moments 

(21) m k {t) with k = j + %l, j, 1 = 0,1... 

Since the total mass is conserved, mo(t) = mo = const; we shall enumerate the rest 
of the moments (1211) by a single index k n , n = 1, 2 . . ., in the increasing order, and 
introduce the notation 

J = {k n : n = 1, 2 . . . } 
for the index set. Also, introduce the normalized moments 

(22) *W = i^- 

where the constant b > will be chosen below depending on a in (J7|). For 6=1 
and k nonnegative integer we have z k (t) = m k (t)/k\ which is the normalization 
appearing in (fl9l) . Also, as is easy to verify by Stirling's formula, 

(23) r(Jfe + &)~Jfe b - 1 r(A; + l), k^oo, 
so the particular choice of b is irrelevant for ffl9|) . 

Given k = k n G J we set z^ k \t) = (^(i), . . . , £fc„_i(0)j a vector with n — 1 
components. 

By the assumptions on mfc(O), we have 

(24) z k (0)<C q^ kEJ, 

for certain constants Co, qo- We shall show that the geometric growth of the nor- 
malized moments is preserved uniformly in time, due to the structure of the system 
of differential inequalities satisfied by Zk(t); this will imply ( fl9l) . The key step is the 
following. 

Lemma 5. Let the sequence of nonnegative functions z k € C 1 ([0, oo)), k G J, satisfy 

(25) z' k (t)<-A k zl + ^(t) + B k F k (z {k \t)), keJ, k>K 
and 

(26) z k {t)<C iq k , keJ, k<h, 

where k* > ~, G\ and q\ are positive constants, A k , B k are positive sequences 
satisfying 

(27) _*>c' 1 1 "^, keJ, k>K, 

B k 

and F k are continuous functions of their arguments such that 

(28) F k (z {k) ) < C 2 q k+ ^ , whenever z k <Cq k , k e J, k > k*. 

Assume that the sequence z k (0) satisfies (I24p . Then z k (t) < Cq k , k e J t > 0, 
where C = max{C , Ci} and q = max{g , <?i}- 
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Proof. Without loss of generality we can assume that G\ = Co and q\ = go- The 
proof will be achieved by induction on k e J, k > k*. For k = k* conditions ([2"oT) 
and (128]) imply 

By a comparison argument for Bernoulli- type ordinary differential equations (cf. [6]), 

(29) z k (t)<max{z k (0),z* k }, 
where z* k is determined from the equation 

A k (z* k ) 1+ £=B k CU k i +§ 
Using condition ff27|) it is easy to verify that z* k < Ci gf, which in view of ( 1291) 
and ([21]) implies z k (t) < C\q\, k = k*. This provides the basis for the induction. 
The induction step follows by repeating the same reasoning for any k > k*. The 
proof is complete. □ 

Next, we shall verify the conditions of Lemma [5] for the sequence of the moments 
corresponding to a solution of the Boltzmann equation. The proof of the time- 
regularity of the moments is standard; we refer the reader to Appendix B for the 
details. We can also use the known property that the moments of every order are 
uniformly bounded in time (part (iii) of Theorem [2]) to deduce (1261) . The main 
difficulty is then to obtain the system (1251) and to make sure that the necessary 
estimates hold for the constants. 

Let us first make some general comments about the time-evolution of the quan- 
tities L d f{v , t) ty(\v | 2 ) dv, where ^ : K.+ — ► R is a convex function. Multiplying 
equation ([2]) by ty(\v | 2 ) and integrating with respect to v we obtain, after standard 
changes of variables, 

(30) -/ f(v,t)V(\v\ 2 )dv = [ [ f(v,t)f(y*,t)Rv(v,v*)dvdv*, 
where 

(v, v*) = \v — vjf (G* (v, f*) - (v, v*)), 

GAv,v*) = \ [ ^(\v:\ 2 ) + ^(\vf))h(^^)da, 
j s ^ 

v'x, v' are defined by (j3J), and 

L*(v,«o = ^(*(H 2 ) + *(KI 2 ))- 

Since the expression for Gy(v,v*) is clearly the most complicated part of (1301) we 
look for a simpler upper bound. This is achieved by means of the following estimate. 

Lemma 6. Let \1/ : K + — > K. be convex and assume that the function h(z) = \{h(z) + 
h{—z)) is nondecreasing on (0,1). Then 

G 9 (V, V,) < LU d _ 2 J V((\V\ 2 + |^| 2 ) li-^ ) h( Z ) (1 - Z 2 )^ dz. 

Proof. See [7, Lemma 1] for the case d = 3; the extension to general n is straight- 
forward. □ 
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Recall that the mass mo and the energy mi are constant for the solution f(v,t). 
We will also use a lower bound for the moments of order a < 1. 

Lemma 7 (Cf. [6] for the case a = 1). The solution of (J2]) satisfies 

f(v*,t)\v-v*\ a dv*>c a fo{v*) \v — tf*| a G?t>*, v G M d , 



/or any a G (0, 1] . 

Proof. By translating the solution f(v*,t) in the velocity space, we can reduce the 
proof to the case v — 0. We will establish the estimates 

(31) m Q (t) > (0), 

for < a < 1. Notice that ^/(^) = — z a is a convex function. Then, by the previous 
computation, and using Lemma [6], 

<{t)> [ [ f(v,t)f(v.,t) \v-v*f(^(\v\ 2 +\v4Y-U\ v \ 2a+ \ v *\ 2a )) dvdv * 

where a a = 2 §\{^Y-) a b{z) (1 — z 2 )^ dz > 1. We shall estimate the integrand 
above in order to obtain an expression involving m a (t) and similar quantities. For 
this we notice that since (x + y) 13 < x 13 + y 13 , for (3 G [0, 1], then 

\v-vjf < (\v\ + < \vf + \v m f. 

Also, 

\v-v»\ p > I \vf - \v»f\ and (\v\ 2 + \v*\ 2 ) a > \ \v\ 2a - \v*\ 2a \. 



Therefore 



\v-v.f (\v\ 2 + \v*\Y-l(M 2a + \v*n) 



> y - \v*f) (\v\ 2a - \v*\ 2a ) ~l(\vf+ \v*f)(\v\ 2a + H 2a ) 

= {\v\ P+2a + \v*f +2a ) - ^±1 (\v\ p \v m \*> + \v\ 2a \v^) 

and we obtain 

m' a (t) > (a a - 1) m m (t) - (a a + 1) mp (t) m a (t) . 

^ 2 2 

In the particular case (3 = 1 we have 

^i(^) > (ai — 1) moffii — (ai + l)m|(t), 

(mo and mi are constants, by the conservation of mass and energy). Therefore, 

f /«i - 1 \n r /ai-ixh 

mi (t) > min< mi (0), — mo mi > > min< 1, — >mi(0), 

I 2 V«| + l /J 1 + J 2 
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since momi > mi (0) 2 . (This is the argument of Bobylev.) To achieve the proof for 

(3 < 1 we iterate this argument, applying it with a = 4^, j = 1 . . ., until > 1. 

Consider first the case of the terminal j, when 

„ ^ * i < a + 1)(3 

= T - • 

In that case 



> {a ao - 1) m m x - (a ao + 1) m m ao (t) 



Therefore, 



i 



m ao W > mm<j m ao (0), ( 7T m o m ] " 

i 



Further, take a x = a — | > 0. Then 

m' ai (t) > (a ai - 1) m m ao (t) - (a ai + 1) m 2Ql m a j 2ai (t), 



i 2_ 

2ai T 2a, 



SO 



m ai (i) > min<! m Ql (0), ( ( — * ) m^ 1 m ao (t) 



a ai + 1 



f- (0) -((^ 



> min^ m Ql (0), ( ( =2L_^ ) ( ) ^'^'m^O) " ^ 



> — r — r m ai (0). 

\a ai + ly \a ao + lJ 

The rest of the proof can be obtained by induction. This establishes fl3Tl) for all 
a G (0, 1] and completes the proof. □ 

In the particular case *&(z) = z k , k > 1, we obtain the following inequalities 

(32) m' fc (i) < / / f(v,t)f(v*,t)R k (v,v*)dvdv*, 

J'R d J'R d 

where 

(33) R k (v, v.) = \\v- v^(a k (\v\ 2 + \v*\ 2 ) k - \v\ 2k - \v*\ 2k ), 
and the constant a k is defined by 

(34) a k = uj d ^ 2 J (^±iyh{ z ){l-z 2 ) d -^dz, 
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satisfies a\ = 1, a k < 1 for k > 1 and is strictly decreasing with increasing fc. We 

1/3 



2 ^kP) 



notice the inequalities |t> — t>*| /3 < + |t>*| 

(M 2 + M 2 ) fc - M 2fe - H 2fc < Q) (kr>*l 2M + M 

where [•] denotes the integer part (cf. [7]). Also, by Lemma El 

/(t>*,t) \v - v*\P dv* > op /o(w*) \v - v*f dv* > u (1 + \vf), 



where i>q is a constant depending on (3 and /o. Using these inequalities in ( 1321) . ( 133 
we obtain 



(35) m' k (t) < -(1 - a k ) v m k+ e(t) + a k S k (t) 

where 



Skit) = ( • ) ( m J+ i (*) m fc-j(*) + m fc-j+i (*) m j(t))- 

3=1 

The crucial estimate for the sum Sfc(t) is provided by the following Lemma. 
Lemma 8. For b > /ixed sei z k (t) = m k (t)/T(k + b), k > 1. T/ien 



5 fc (t)<C 6 r(A; + | + 26)Z fc (t), fc > 1, 



where 



(36) Z fc (t) = max {z j+ p (t) z k -j{t), Zj{t) z k _ j+ M)} 
and Cb is a constant depending on b. 

Proof. See [7, Lemma 4]. □ 

Proof of Theorem^ Using the interpolation inequality m k+ p(t) > m 2k m k (t) 1+ 2k 
and Lemma [8] we derive from (|35|) the inequalities 

(37) z' k (t) < -(1 - a k ) u mo - T(k + b)& z\ +2 (t) + a k C b E^ + l^ Zfc (t). 

Notice that for k G J, k > 1 + | the term Zfc(t) is of the form F k (z( k \t)) as 
in Lemma since the highest order of moment entering (1361) is fc — 1 + | . It is 
also clear the the function F k defined in this way is a continuous function of its 
arguments. Thus, we can identify (|37|) with (125]) by setting 

s- §_ „ „ r(fc + f + 26) 



(38) A k = (l-a k )u m 2k T{k + b)^ } B k = a k Q 



T(k + b) 



We would like to apply Lemma [5] to the sequence of functions z k (t). It remains 
to verify that the sequences of constants A k and B k appearing in (J3"7|) satisfy (j2j 
To this end we show that 

(39) 4^ > c o> fc > 
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for any k* > 1 + | and a sufficiently small Co; then (127)) would follow by choosing 
Co = Ci small enough and go = Qi large enough in ( 124"1) . ( 126)) . Indeed, using ( 1231 ). 



r(fc + f + 2b) 



oo. 



(40) r(Jfe + b)~™ ~ A;"/ 2 and r ( fc 2 +& ) ~ fc > fc 

To estimate the constant in (1571) we recall that b(z) < C (1 — £ 2 )~ Q , a < cf — 1 
and setting in (154")) . s = ^rp, e = n — 1 — a>0, we have 

• i 



(41) 



a fc = C2- 1+£ / s fe - 1+ i(I - s)- 1+ i = C2- 1+£ B{k + £ £ 



o 



2' 2' 



r(A; + £) 



We fix < b < e/2; the corresponding constants satisfy the inequalities 

(42) A fc >AA;i, B k <Bk2 +b -^, k > h m , 

where fc* > 1 + f , and A and 5 are absolute constants which can be estimated based 



on ([38]) and the asymptotic relations ([40]) and (j4lj) . From ([42]) we obtain ([59]) by 
choosing c = AB -1 ki . 

We conclude the proof of Theorem [4] by applying Lemma [5j □ 

3. Comparison principle for the Boltzmann equation 

In this section we discuss the important technique of comparison that will allow us 
to obtain pointwise estimates of the solutions. The crucial property of the Boltzmann 
equation used here is a certain monotonicity of a linear Boltzmann semigroup. The 
argument is roughly as follows: if / is a solution of (CQ), f\t=o = fo, and g is 
sufficiently regular and satisfies 

(43) (d t + vV x )g>Q(f,g), g\ t=0 = g , 
then u = f — g is a solution of 

(44) (d t + vV x )u<Q(f,u), u\ t=0 = u , 

where Uq = fo — g^. We will show that if / is nonnegative (and satisfies certain 
minimal regularity conditions), then solutions of (jUj) satisfy the order-preserving 
property, 

(45) if m < then u < 

(zero on the right-hand side can be replaced by any other solution u of ( 144)) ). This 
translates into the following estimate (comparison principle): 

(46) if fo < go and g satisfies ( |4"51) . then / < g. 

By reversing all inequalities we obtain a similar comparison principle that yields 
lower bounds of solutions. 

Of course, the above scheme has to be implemented with suitable modifications. 
For instance, since apriori only limited information about / is available we will 
require that g satisfies ( 1151) for a class of functions / (defined by the available apriori 
estimates). Another important refinement is to apply the estimate (j46p locally (in 
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the case of Theorem [3l to a "high- velocity tail" {\v\ > R}) since global bounds in 
all of the (v , t)-space cannot be generally obtained by this technique. We refer to 
Proposition [9] and the proof of Theorem [3] given below for the necessary details. In 
Theorem [TTJ we will give a rigorous statement of (H6|) in application to a general 
class of weak solutions of ([TJ) in the sense of DiPerna and Lions [14, 28] . 

The basic approach leading to applications of (]4"6l) originated in the work by 
one of the authors [34, Sec. 6.2] in the context of lower bounds for the spatially- 
homogeneous equation without angular cutoff. It was also used to obtain lower 
bounds for solutions in a model describing inelastic collisions [18]. Compared to 
these earlier versions we do not require in (j4"6l) any differentiability in the f-variable, 
and we make more precise the minimal regularity conditions on /. It is interesting to 
compare the present technique with other methods based on monotonicity applied 
to the Boltzmann equation, in particular the one by Kaniel and Shinbrot [25] (see 
also [21,24]) and the pointwise estimates by Vedenjapin [33] (the result in the latter 
paper follows from our approach using g = e C( - 1+ *- ) ). The monotonicity property ex- 
pressed by ( l4"5l) has also an important relation to the concept of dissipative solutions 
introduced by P.-L. Lions [28]. 

We first explain the way to obtain (T4"5]) . The bilinear form in ([4*3]) . (|44"I) is defined 

by 

(47) Q{f,u){x,v,t) = / {flu' - f*u)B{v - v*,a) do dv*, 

Jm. d J s 6 -- 1 

where as usual, fl = f{x, v' # , t),u' — u{x, v', t), /* = f{x, t>*, t), u = u{x, v, t). At this 
point we do not need to assume the kernel B to satisfy (j5j)-([7j); the argument goes 
through for a more general class of kernels with the usual symmetries, as described 
in [14], for instance. 

To illustrate the general principle, consider first the case of equality in (1441) . Given 
T > we fix the function / : T x R d x [0, T] — > M + , which we assume to be smooth 
in {x,t), bounded and rapidly decaying for \v\ large. We also assume that for every 
Uo G D C L l {T d x W 1 ) the initial- value problem 

(48) (dt + v- V x )u = Q{f,u), u\ t=0 = u , 

has a unique solution u G C([0,T]; L x {T d x M d )), with enough regularity so that 

(49) Q + {f, \u\), Q~{f, \u\) E L\T d x R d x [0, T}). 
Thus, we have a well-defined flow map (or a semigroup) 

$t : D 3 u i-> u(t, v) e ^(T d x R d ), t E [0,T]. 

The map <3?t can be seen to satisfy the following nonexpansive property: for any 
w , u G D, 

(50) / / \§ t {uo) - $t(«o)| dvdx < / / \u Q — u \dvdx, te[0,T]. 

Indeed, set w = $t(wo) — &t{uo); then 

{d t + v-V x )w = Q{f,w) on T d xE d x(0,T) 
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in the sense of distributions, and Q(f,w) G L 1 by our assumptions. By a standard 
argument, Vt G [0,T], for a. a. (a;,t>) the function : s h> w(x — (t — s)v,v,s), 
s G [0,T], is absolutely continuous, and we can apply the chain rule (see Appendix 
A) to obtain 

(51) — = Q{f,wf signw tt , s G (0,T), 

where Q(f,w)$ is defined similarly to wK Integrating with respect to s G (0, t) and 
(x, v ) G T d x ]R d we obtain, after standard changes of variables, 

/ / \w(x, v, t)\ dvdx — I I \wo\dvdx+ / / / Q(f,w) sign wdvdxds 

JR d JR d jR d jR d JO jR d JR d 

where Wq = Uq — u . We further notice that the bilinear collision term (|47p satisfies 

(52) / Q{f,u) signudv < 0, 

JR d 

for every / > and every u so that Q + (f, \u\), Q~(f, \u\) G L 1 . This follows 
immediately from the weak form 

Q(f, u) sign u dv = / / / f*u (sign u' — sign u) B da dv*dv 

JR d JR d J S 1 *- 1 

by noticing that w (signw' — signw) < 0. 

The same approach can be followed to obtain (145j) . Indeed, we have by (1521) and 
the mass conservation 

(53) 



/ Q{f,u)±{signu + l)dv <0, 

JR d 



where |(signw + 1) is the a. e. derivative of the Lipschitz-continuous function u + 
m&x{u, 0}. We then have 

± u { = Q(f,uf l(signu + lf, se(0,T), 



and the integration yields 



/ / u 0+ dvdx, ie[0, T], 

JR d JR d 



u + (x, v, t) dv dx < 
which implies (1451) for 

Remark. Relation (145]) can be restated as the order-preserving property of <i\: 

(54) Vu , n G £>, u < «o implies $ t (« Q ) < $t(« ), * e [0,T]. 

In fact, the equivalence of ( 1541) and ( |50|) follows from a general principle applied to 
(nonlinear) maps that preserve integral, as described by Crandall and Tartar [11]. 
Inequality (145!) (or (|54l) ) can then be seen as a consequence of the results in [11], 
the preservation of the mass f Rd L d / along solutions of (jHJ) , and (150]) . 

The following localized version of the order-preserving property will be useful for 
the comparison argument. 
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Proposition 9. Let f,u e C([0, T\\ L l (Y d x R d )) satisfy 

/>0; d t u + v -V x u, Q + (f,u), Q~(f,u) G L 1 ; u\ t=0 = u <0, 
and assume that for a certain (measurable) set U C T d x IR d x (0, T) 7 

c^n + n • V ' x u — Q(f, u) < on Z7, 

and 

u < on f/ c := (T d x R d x (0, T)) \ U. 
Then u(t, ■, •) < a. e. on T d x ]R d , /or every t G [0,T]. 

Proof. Let -D(n) = 9jM + f ■ V x n. We obtain by arguing as above, 

u + (x,v,t) dv dx — I I u + (x,v,0) dv dx 



/ / / D(u) | (sign u + 1) dv ds. 



< 




We have u + \ t=0 = 0; also |(signn + 1) = whenever u < and -D(n) = outside of 
a set of zero measure in {u = 0}. Therefore, setting U t = {(x, v , s) G U : s < t} we 
have 

u + {x,v,t) dv dx = III D{u)dxdvds 
J J J u t 

Q(f,u)dxdvds = Q{fi u ) |(signn + 1) dxdv ds < 0, 

t Jo Jtt d Jtt d 

for every t G [0,T], where we used the dissipative property fl53l) . This shows that 
u(t, •) < almost everywhere. □ 

Proposition [9] is sufficient to formulate the comparison principle in the generality 
required for Theorem [TJ We will, however, give a more general statement that 
applies to weak solutions in the spatially inhomogeneous case. In the definition of 
weak solutions one has to account for the fact that the bound 

Q(f)eLl c (T d xR d x(0,+oo)) 

is generally not available, and one has to define solutions in a sense that is weaker 
than distributional. The simplest way to state the definition is to require that / > 0, 
/ G C([0,T]; Ll v ), Q ± {f)/{1 +/) G L\ oc and the renormalized form 

(d t + v-V x ) log(l + /)=W)/(l + /) 
holds in the sense of distributions, cf. [14]. Such solutions are known as renormalized. 
This concept can be further refined as follows, cf. [28]. 

Definition 10. We say that a renormalized solution / is dissipative if f\v | 2 G 
L°°([0, T]; L x xv ) and for every sufficiently regular function g : T d x R d x [0, T] -> R, 

(55) <9 t / \f-g\dv + dw x f \f-g\vdv< [ (Q(f,g)-D(g))sign(f-g)dv, 

JR d JR d JR d 

in the sense of distributions, where D(g) = (d t + v • V x )g, and sign(0) is assigned an 
arbitrary value in [—1, 1]. 
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Remark. In the above definition "sufficiently regular" precisely means that g G 
C([0,T];Ll v ), g\v\ 2 G D(g) G L\ vt and that for any / G C{%T\-L\ V ) 

such that f\v | 2 G L^°(L^), Q + (f, \g\), Q (f, \g\) G L l xvt (these conditions can be 
made more explicit, see [28] for details). 

The formal motivation for the definition of dissipative solutions is clear: the right- 
hand side of the Boltzmann equation can be written as 

Q(f) = Q(fJ-g) + Q(f,g): 

so we have 

(d t + v- V x )(f -g) = Q(f, f-g) + Q(f, g) - D(g). 
Multiplying the above equation by sign(/ — g) and using relation fl52l) (note that 
/ > 0) we see that every sufficiently regular solution of §1§ should satisfy fl55|) . 

Dissipative solutions are known to exist globally in time, for a quite general class of 
initial data. In fact, in [28] Lions established a large class of "dissipation inequalities" 
similar to (1531) that hold for renormalized solutions of ([IT) . Such solutions can also 
be constructed so that the local mass conservation law, 

(56) d t [ fdv + div x [ fvdv = 0, 

jR d JR d 

holds in the sense of distributions. However they need not generally satisfy the 
conditions Q + (f), Q~(f) G L[ oc . 

Using the order-preserving property of Proposition [9] we establish the following 
comparison principle for dissipative solutions of the nonlinear Boltzmann equation. 

Theorem 11. Let f G C([0, T]; L l (Y d x R d )) be a dissipative solution of (pQ) and 
let g be a sufficiently regular function, such that f\t=o < g\t=o, 

d t g + v-V x g-Q{f,g) > on U 

and f < g on U c , where U is a measurable subset ofT d x M. d x [0,T]. Then f < g 
almost everywhere on T d x IR d , for every t G [0, T\. 

Remark. It is natural to call g a (localized) upper barrier. By reversing all inequal- 
ities in the above formulation one can also obtain a similar comparison principle for 
the lower barrier. 

Proof. We use the notation D(g) = d t g + v ■ V x g, so that 

d t gdv + div x gvdv= / D(g)dv, 

JR d J~R d JR d 

in the sense of distributions. Using the mass conservation (1561) and the identity 

(f-g) + = l(\f-g\ + (f-g)) 

we obtain, by combining the above relations with ( |55l) . 
dt (/ - g)+dv + div x / (f-g) + vdv 

JW. d JR d 

<\ I {Q{f,g)-D{g))siga(f-g)dv-l [ D(g)dv. 

/ JR d A JR d 



MAXWELLIAN BOUNDS FOR THE BOLTZMANN EQUATION 17 

Since Q ± (f, \g\) are integrable, we have L d Q(f, g) dv = 0, a. e. (x, t), and therefore, 
9t (f - g)+dv + div x / (f-g)+vdv 



(57) 



< 



/ (Q(/^)-D(^))i(sign(/-(/) + l)dT;. 



We can choose sign(0) = — 1 in ( 1571) to avoid estimating the integral over the set 
{/ = g}. Since (/ — g) + v G L 1 {T d x R d x [0,T]) we can integrate over x and t to 
obtain 



(5S 



(/ - g)+(x,v,t)dvdx < / (f - g) + (x,v,0) dvdx 

+ II lu (W^)-D(g))dxdvds<0, 

where £7 t = {(x, v, s) E U : s < t} and we used that |(sign(/ — g) + 1) vanishes for 
f < g and that — -D((?) < on f/ t . The inequality in fl58l) implies that f < g, 

a. e. (ar, u) G T d x M ', for every t G [0, T\. □ 

Theorem [TTl is a crucial ingredient in the proof of Theorem [3l which we give below. 

Proof of Theorem^ To apply Theorem [TTl we set U = {(x,v,t) : |u| > R}, where 
R will be chosen large enough, and g(x,v,t) = M(v), where M(v) = e~ a ^ +c , 
< a < ai is fixed and c > Co will be chosen sufficiently large, depending on R. To 
prove that g can be used as a barrier for the solution on U we need to verify the 
inequality 

(59) Q + (f,g) (a;, v,t) < Q-(f,g)(x,v,t), (x,t) G T d x [0,T], \v\ > R. 
First notice that, by elementary inequalities, 

Q~(f,g)(x,v,t) = M(v) / f(x,v*,t) \v -v+fdv* 



>M(v)(p \vf- f(x,v„t)\v,\ p dv< 



where po is the constant in ffT3]) . The last term can be controlled using the estimate 
for the integral of //Mi from (|14p as follows, 

f(x,v*,t) \v*fdv* < L [ ^ ?J V ,*'l dv* < LCx, 

where L = max y 13 e~ aiy2+Cl . Thus, we have 

y>o 

Q- (f, g) (x, v, t) > M(v) { Po \vf -Ld). 

The control of the "gain" term is more technical; we establish below in Lemma 
the estimate 

(60) Q+(/, g)(x, v,t)<C(l + \vf- £ ) M(v), 
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where e = min{/3, n — 1 — a} > 0. This implies that (l59i) holds if we set R to be the 
largest root of the equation 

C + LC 1 + Cy^ - poyf 3 = 0. 

Finally, we take c = aR 2 + logC , where C is the constant in ffl4"I) ; then it is easy 
to verify that 

(61) f(x, v, t)<C < M(v), (x, t) e T d x [0, T], \v\ < R. 

The conditions < a < a x < a and c > c guarantee that we have f(x, v, 0) < 
M(t>). Together with the inequalities fl59|) and flBT]) this allows us to use Theorem [TT1 
to conclude. □ 



4. A WEIGHTED ESTIMATE FOR THE "GAIN" OPERATOR 

To complete the proof of Theorem [3] we prove the following weighted estimate of 
the linear "gain" operator. The main technique is based on Carleman's form of the 
"gain" term (see Appendix C). 



Lemma 12. Let B : R d x S^ 1 -> 

B(u,a) <C(1 + \uf) 



, n > 2, be a measurable function that satisfies 

l{costf>0}, COS?? = TTT, 



sin^l' 



M 



where (3 > and a < n — 1. Define 

Q + (f,g)(v)=[ [ f: g 'B(v-v*,o-)dadv*, 

and set M(v) = e~ a ^ 2 , a > 0; iu e (u) = 1 + |f l' 3 ^ 5 , where e = min{/3, n — 1 — a} > 0. 
Then 



(62) 



Q + (/,M) 



< C 



L°° 



[We 

M 



L 1 



where C is an explicitly computable constant depending on n, a, (3 and a. 

Remark. For B satisfying the estimate with a = (for example, the kernel B for 
hard spheres in three dimensions) we have e = (3 for all /3 < d — 1 and the weight 
w e (v) is constant. The estimate of the Lemma then takes a particularly simple form, 



Q + (f,M) 


< c 


f 




M 




M 


Li 



For the quadratic "gain" term this implies the estimate 



Q + (f) 

M 



< C 



f_ 

M 



Proof. By the Carleman representation formula, 

rd-i I f« 



Q + (f,M)(v) 



f_ 

M 



M(V') B{V - V *> a ) dTT, 



v - V' # \ J E , V- V* 

where E vv > t is the hyperplane 

{v 1 G M. d : (v- v') ■ (v - vi) = 0}, 



n-2 
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and dii v i denotes the usual Lebesgue measure on E vv /^ . We then have 

Q + (f,M)(v) f ffo) v( 

(63) TTTS = / —r, — r K(v, v ) UV , 

y ' M(v) J Rd M«) v * ; 

where 

(64) K(v,vi) = [ M(v.) B{V ~ V Z^ dir & , 

\ v - K\ Je , \ v ~ v *\ 



and we used that, by the energy conservation, 

M(t/)M«) 



MO* 



M(v) 

Note that in (1641) the variables i>* and a are expressed through v, v[ and v' as follows, 



u* = v'+v' — v, n 



v' — v' 



\v' — v'J 

Now to establish the Lemma it suffices to verify the inequality 
(65) K(vX) < C(l + \v-v'f- £ ). 

Indeed, since 

l + \v-vl\P- £ <(l + \vf- £ )(l + \v'f- £ ), 
then (1631) and ( 1651) imply 

Q + (f,M)(v) < C{l + \v\P-*)M(v) [ *(l + |<h*; 

which is equivalent to (l6"2"j) . 

In the remainder of the proof we will therefore verify fl65|) . Using the identity 

(u - v*) • (?/ - u*) = |u - -u^ 2 - |u - i/| 2 

for f ' G E vv i^ and recalling that — f*, a) vanishes for (v — vj) ■ a < we see that 
the integration in AM]) can be restricted to the disk 

D v< = E v< n {V G R d : |u - v'J < \v - v'\}. 

We notice that for v' G Dvv' t , 

l tan 2l = F^' "-2' 
where i? is the angle between the vectors v — v * and a. This implies 

1 1 \v — v'J 



| sin i? I 2 — vj 
Thus, if(i>, i>*) < CK(v, v'J) , where 

( v , <) = / M(^) , 1 , * , - , dn v , 

\v-vJ l ~ a J D , K ' \v-vj n ~ 2 k-vJ a 

To estimate the above expression we consider two cases. 
Case a) \v — v'J < 1. Since for v' G D vv >^ 

\v — v'J < \v — vj < v2 \v — v'J 
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we have 1 + 1 1> — f*| /3 <l + 2^/ 2 and 
Therefore, 



v-vj 2 - n < \v-v'J 2 - n . 



< n — ;i»-i- a / n — 

« - v *\ Jd , K~ v *\ 

Since M(v*) < 1 the last integral is estimated above by 



f — 



if c? — 1 — a > 0, i. e. a < d—1. Here 00^-2 is the measure of the (n — 2)-dimensional 
unit sphere. This implies the estimate 

a — 1 — a 

Case b) \v — v'*\ > 1. Then 

1 + \v - v*f <2\v- vj[ p < 2 1+ 2 \ v - vif, 
and we obtain, similarly to the previous case, 

K(v,v'J < - — — s / M(v,) — -d-K v >. 

Since M(t>*) is a radially decreasing function of v* G M d , and so is |f*|~ a , 

< / |u ; |- Q rfw+ / M(to) dw = Ud ~ 2 +(-)~, 
J\w\<i JR-i- 1 d—l — a \a/ 

where M(w) = e~ a ^ 2 , w G R^ -1 . Since \v — v 'Ji 3 + a - n + 1 < \ v — Vj>t \P-£ this establishes 
the required estimate for Case b). □ 

Appendix A: Some properties of weakly differentiable functions 

Let AC [a, b] denote the class of absolutely continuous real-valued functions defined 
on an interval [a, b}. Given / G AC [a, b] we set [c, d] = /([a, b]) and use the notation 
Lip [c, d] for the class of all Lipschitz continuous functions defined on [c,d\. Every 
function (3 G Lip[c, d] is differentiable (in the classical sense) almost everywhere on 
(c, d); we agree to extend this derivative to a function (3' defined everywhere on [c, d] 
by assigning arbitrary finite values at the points where (3 is not differentiable. The 
function (3' also coincides with the weak derivative of (3 almost everywhere on (c, d) 
The following chain rule was used in the arguments in Section [31 

Proposition 13. Let f G AC [a, b] and (3 G Lip[c, d]. Then /3 o / g AC [a, b] and 

09 o /)' = {(3' o /) /', 

almost everywhere on (a,b). 
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Remark. 1) The seeming ambiguity in the above formulation occuring since (3' o f 
can assume arbitrarily assigned values on a set of positive measure is resolved by 
observing that whenever this happens then /' vanishes, except on a set of measure 
zero (see the proof below). 2) For the purposes of Section [3] we only need the chain 
rule for /3(y) = \y\ and (3(y) = y + ; these cases are covered in [17], and the proof for 
the case of piecewise-C 1 functions (3 can be found in [20]. We include a short proof 
that applies to the general case to make the presentation in Section [3] self-contained. 

Proof. By the definition of absolutely continuous functions, 

Ve > 35 >0 such that VneN, V {(xj, yj) C [a, b] : j = 1, . . . , n}, 
a disjoint family, 

n n 

3=1 i=i 

Clearly then, since 

\P{f{yj))-P{f{xi))\<L\f{y j )-f{x j )\, 

where L is the Lipschitz constant of /3, the composition /3o / is absolutely continuous 
on [a, b}. By Lebesgue's differentiation theorem, / and f3o f are differentiable in the 
classical sense on a set with complement of measure zero in (a,b). Pick x G (a, b) 
from this set. We will consider two cases, depending on whether (3 is differentiable 
at f(x) or not. In the first case we have 

™ f(x + h)-f(x) ™ ft 

Let us further take A to be the set of y such that (3 is not differentiable at f(y). 
We claim that f'(x) vanishes for x G A, except perhaps on a set of zero Lebesgue 
measure. Indeed, let B = {y G A : > 0}; then 



B= UB n , B n = {yeB: \f(z) - f(y)\ > ^ for \z - y\ < J }. 



oo 

u 

n=l 

We prove the claim by showing that every set B n has zero measure. 

Fix an n G N. Since (3 is Lipschitz, we know that f(A) is a set of measure zero. 
Given e > we can then choose the intervals Ij, j = 1, . . ., such that 



f{A) C U^/j and 2^ l^-j < e. 



Let J be an interval of length and let D = B n D J, -Dj = HD. Then, from 

the definition of \Dj\ < n\Ij\; therefore, \D\ < ne and \B n \ < n 2 \b — a\e. Since 
e is arbitrary this shows that \B n \ = 0. 

We now have that for a. a. x G A 

+ /?(/(*)) <l /(* + /0 - /(x) 



h 



h 
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for \h\ small enough, so (/3o/)'(x) = and fl'(f(x))f'(x) = 0. This proves the claim 
of the Lemma for a. a. x G (a, 6). □ 



Appendix B: Time regularity for the spatially homogeneous 

BOLTZMANN EQUATION 

We show that the solution of the Boltzmann equation ([2]) under the conditions of 
Theorem [Tj is smooth with respect to time, together with its moments of any order. 

For k > we introduce the following weighted Lebesgue spaces 
(66) Ll(R d ) = \f e L^W 1 ) : f \f\ (1 + \v\ 2 ) k dv < +oo) 

L JR d J 

with the norms defined by the integrals appearing in fl66|) . The regularity result that 
we used in Section [2] is the following. 

Proposition 14. Let f be the unique solution of the Boltzmann equation ([2]) that 
preserves the total mass and energy. Assume that fo e Ll(R. d ), k > 1 + f . Then 
f e C 1 ([0,+oo);Lj(R d )) for anyp<k-\. 

The proof of Proposition [TD depends on the following continuity property of the 
nonlinear operator Q{f). 

Lemma 15. Let the pair of positive numbers (k,p) satisfy k > p+ |. Then Q(f) 
is continuous on Ll(R d ) as a mapping L£(R ) — > Lp(M. d ). Moreover, we have the 
following Holder estimate for any f,gE L\(R d ) 

II W) - Q{g)\\ L i < c p (\\f - gW 1 ^ + ||/ - ^|| L i) , 

where the constant C p depends on p and on the upper bound of the L\-norms of f 
and g. 

Proof. Using the weak form of Q(f) and Q(g) we compute 

/ \Q(f)-Q(g)\^ + \v\ 2 ) p dv 

= f f f (ff*-gg*)B(v-v^a)(sign(Q(fy-Q( g y)(l + \vfr 

JR d JR d JS"- 1 V 

- sign(Q(/) - Q(g)) (1 + \v\ 2 ) p ^j da dv dv, 

<2 P+1 [ [ \ff, - ggJ \v - v*f ((1 + \v\ 2 ) p + (1 + \vJ 2 ) p ) dvdv* 
Jm d Jm d 

Since 

\v - v*f (i + \v\ 2 y < (i + m 2 ) f (i + \v\ 2 y + (i + \v\ 2 y + -i 

<2({l + \v\ 2 Y + *> +{i + \vA 2 y + -*) 
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and |//, -gg*\ <\\f-g\ \ f* + 9*\ + \ \f +g\ \f* - g*\, we obtain 

\\Q(f) - Q(g)hi 

<2 p+3 [ [ \f + g\\f*-g*\ ((l + \v\ 2 ) p+ % + (l + \v*\ 2 ) p+ %)dvdv* 
<2 p+3 \\f + gWa (\\f - g\\ L i B +\\f-g\\v). 



We use the interpolation inequality (1201) with k\ — p + | to get 



\\f-g\W <\\f-g\\ L , k \\f-g\\A 



, P+% E±3L 
1 '- " - " k 



Substituting this bound into the previous estimate we obtain the Holder estimate 
stated in the Lemma. This completes the proof. □ 



Proof of Proposition 14 ■ We fix T > 0. By the results of Arkeryd and Elmroth [1,15] 



(see part (iii) of Theorem [2]), / belongs to L°°([0, +oo); L\{R d )). By Lemma [T5| 
(67) (l + \v\ 2 YQ(f)eL 1 ((0,T)xR d ), for p < k - | 



The mild form of (J2J), together with the regularity condition fIBTj) imply that / is 
weakly differentiable and dtf = Q{f) in the sense of distributions on (0, T) x M. d . 
Hence, 

and therefore (cf. [16, p. 286]), / £ C([0, T]; L^(R d )) . By the continuity of Q(f) 
established in Lemma [151 it follows that d t f £ C([0, T]; Lp(IR d )), where <9 f / is the 
weak time-derivative of /. It is then easy to verify directly that / is strongly 
differentiable on (0,T) with values in Lp(IR d ), and its derivative is continuous on 
[0, T\. Since T is arbitrary, we obtain the conclusion of the Lemma. □ 

Remark. As a consequence of Proposition if the moments of all orders are finite 
initially then they are continuously differentiable functions of time. By iterating 
the argument we used in the proof above one can show that in fact then / £ 
C°°([0, oo); Ll(R d )), for any k > 0. 

Appendix C: Carleman's representation 

Lemma 16. Let Q + (f,g) be defined by (Q and let f = f(v) and g = g(v), v £ R d 
be smooth functions, decaying rapidly at infinity. Then 



/K) r g( l /)B(2v-v'-v'„^ [ : 



Q+(/, g)(v) = 2" / 7 - _ * Ar„, 44, 

JM d l w v *\ JE v , \ v v *\ 

where E v y t is the hyperplane {v' £ lR d | (v' — v) ■ (i>£ — f ) = 0} and c?7iv denotes the 
Lebesgue measure on this hyperplane. 
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Proof. Using the change of variables u = v — v*, and recalling the definition of the 
delta function of a quadratic form, see [19], we have 

(68) Q + (f,9)(v)= [ [ f(vi)g(v')B(u,k)6(^)dkdu, 



where v' = v — u + | (u + \u\k) and v'^ = v — | (u + \u\k). We further set z = 
— \{u + \u\k); for every u fixed this defines a linear map k i— > z with determinant 

•M 



\k\ 


2 -l 


+ M 


2 




It 


2 2z- 


(z + u) 




2 


2| 


u\ 


2 






u\ 


2 



(^■) . We also have 

2z + u 

k = : — : — and 

\u\ 

With this change of variables the integral in (|68j) can be written as 

/ / (^) d f(v + z)g(v-u-z)B(u,-^)5(^l)dzdu. 

JW JM. d ill 

We set y = —z — u; then \u\ = \y + z\ and 8 ( 2z u ^ ) = ^ y+ ^ S(z ■ y). Further, for 
any test function (p, 

/ S(z ■ y) ip(y) dy = \z\~ 1 j ip(y) d-!r y , 

JR d Jzy=0 

where dn y is the Lebesgue measure on the hyperplane {y : z ■ y = 0}. This yields 

Q + (f,g)(v) 

= [ [ f(v + z)g(v + y) \z\~ l \y + z\ n ~ 2 B(-y - z, dn y dz 

Jz£R d Jy-z=0 

We now return to the original notations v'^ = v + z, v' = v + y and perform the 
corresponding changes of variables to obtain the expression for Q + (f,g) stated in 
the Lemma. □ 

Remark. The above result takes a particularly simple form in the case of the hard- 
sphere model in M 3 ; in that case B(v — i>*, a) = ^-\v — v*\ and 

Q+(f,g)(v)= [ -f^L f g(v')d7T v/ dvi. 
J R 3 7T\V-V'\ J e 
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